Multistability of free 
spontaneously-curved anisotropic strips 

By L. GiOMi AND L. Mahadevan 

School of Engineering and Applied Sciences, and Department of Physics, Harvard 
University, Pierce Hall 29 Oxford Street Cambridge, MA 02138, USA. 

Multistable structures are objects with more than one stable conformation, exem- 
phfied by the simple switch. Continuum versions are often elastic composite plates 
or shells, such as the common measuring tape or the slap bracelet, both of which 
exhibit two stable configurations: rolled and unrolled. Here we consider the en- 
ergy landscape of a general class of multistable anisotropic strips with spontaneous 
Gaussian curvature. We show that while strips with non-zero Gaussian curvature 
can be bistable, strips with positive spontaneous curvature are always bistable, 
independent of the elastic moduli, strips of spontaneous negative curvature are 
bistable only in the presence of spontaneous twist and when certain conditions on 
the relative stiffness of the strip in tension and shear are satisfied. Furthermore, 
anisotropic strips can become tristable when their bending rigidity is small. Our 
study complements and extends the theory of multistability in anisotropic shells 
and suggests new design criteria for these structures. 
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1. Introduction 

The notion of elastic multistability has drawn considerable attention in the past 
few years owing to the potential for the design of smart structures. A multistable 
structure is an elastic object (typically a fiber- reinforced composite) that exhibits 
more than one equilibrium conformation and can thus be arranged in a variety of 
shapes without inducing permanent deformations and with no need of mechanical 
hinges. Everyday examples of these include snapping hair-clips, the slap-bracelet, 
and various jumping toys. Unlike more conventional engineering structures, where 
large deformations must be accompanied by large forces, multistable objects can 
switch between shapes using a small actuation force. This latter feature has made 
multistable structures promising candidates for the realization of a new generation 
of adaptive devices, in which these "morphing" capabilities, combined with limited 
actuation, allow for switchable, controllable conformational changes. Examples of 
this new generation of devices are shape-changing mirrors for adaptive focusing in 
optical systems, morphing aircraft structures that can continuously readjust their 
shape to optimize aerodynamic function (Abdulrahim et al. 2005, Mattioni et al. 
2007). 

The quest for a theoretical description of elastic multistability dates back to the 
1920 when deployable devices, such as the measuring tape, first made their appear- 
ance. In this context, it is amusing to note that A. E. H. Love, the author of the 
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classical treatise on elasticity (Love 1927), was once challenged to explain why the 
longitudinal curvature of a spring-steel measuring tape appears to match exactly 
its transverse curvature (Petroski 2004). He did not; indeed such an explanation 
was offered only decades later by Rimrott (1965) and reviewed and distilled into 
an elegant solution by Calladine (1988). Over the last decade the work of Iqbal et 
al (2000), Galletly & Guest (2004a, 2004b) and Guest & Pellegrino (2006) simpli- 
fied and clarified the elasticity of bistable objects, unraveling the interplay between 
anisotropy and spontaneous curvature in determining the shape of bistable plates. 
Guest & Pellegrino (2006) used a beam model of a strip of infinite length a finite 
width endowed with constant spontaneous curvature along the transverse direc- 
tion. Assuming inextensibility, these authors showed that orthotropic strips exhibit 
a secondary equilibrium conformation in addition to the base configuration. The 
stability of such a secondary minimum of the elastic energy depends on the relative 
magnitude of the bending stiffness in the longitudinal and transverse directions as 
well as the twisting stiffness of the strip. Later, Seffen (2007) considered an ex- 
tensible elliptical plate with free boundaries and non-zero, constant spontaneous 
principal curvatures and showed that even isotropic shells might be bistable. More 
recently, Vidoli & Maurini (2009) extended Seffen's uniform curvature model and 
showed that two-dimensional orthotropic plates with initial shallow double curva- 
ture are in fact instable in some range of elastic moduli and spontaneous curvatures, 
while a one parameter family of continuously variable neutrally stable shapes was 
reported by Seffen & Guest (2011) in the context of prestressed shells. 

In this paper we complement these different results to understand the phase 
space for multistability in spontaneously curved elastic strips, but lift the assump- 
tions of inextensibility and uniform curvature, with the goal of analyzing an entire 
range of morphing scenarios (including tristability) without compromising the sim- 
plicity of the analytical treatment. Our system consists of a free anisotropic elastic 
strip with spontaneous double curvature and, in general, spontaneous twist. In this 
setting, we look for a general solution without assuming deformations to be inex- 
tensible or the curvature to be constant across the strip. This latter feature, in 
particular, is what makes tristability possible even in the reduced dimensionality 
of strip-like plates and allow us to map a "phase-diagram" for the existence of bl- 
and tristability. In section |2] we present the fundamental equations describing gen- 
eral anisotropic elastic strips following the approach of Mansfield (1973), Reissner 
(1992) and most recently, Galletly & Guest (2004b), who considered the specific 
case of composite bistable tubes. In section [3] we specialize our analysis to the 
case of orthotropic strips and discuss the associated energy landscape. In section 
[4] we briefiy discuss the case of strips with coupling between stretching and bend- 
ing, highlighting the differences and similarities with orthotropic strips. Section [5] 
concludes the paper with a focus on open problems and applications. 



2. Elasticity of anisotropic strips 

(a) Formulation 

We start by considering a long rectangular plate whose length is much larger 
than its width 2a, with the x the coordinate along the direction parallel to the short 
edge and y the coordinate along the direction parallel to the long edge (Figure fl]). 
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Figure 1. Schematic representation of a two-dimensional strip of thickness t and width 2a. 
The transverse and longitudinal directions are labelled as x and y respectively. For strips 
that are six times longer than they are wider, it is sufficient to consider them as effectively 
infinite along the y direction. 

In practice we can consider the strip infinite along the y direction to avoid dealing 
with end conditions. As noted by Mansifeld (1973), strips with aspect ratio five or 
more are already well described by this approximation. 

The general constitutive equations for anisotropic plates can be expressed in 
terms of the ABD matrix (Ashton & Whitney 1970; Jones 1999), that relates the in- 
plane stresses N = {N^^ Ny, N^y) and out-of-plane moments M = {M^, My, M^y) 
with the strain e = {cx^ey, Cxy) and the difference % = — c between the middle- 
surface curvature^ = {kxi i^y, i^xy) and their preferred (spontaneous) values c = 

i^^x •) ^y •) ^xy^ 1 SO that I 



N 
M 



D 



(2.1) 



where A represents the extensional stiffness matrix, D is the bending stiffness 
matrix and B embodies the possible coupling between stretching and bending due 
to material anisotropy. 

In terms of the strain, curvature and the material properties of the shell, the 
elastic energy density of a general anisotropic plate is given by (Mansfield 1989): 

1.. ^ . ^2.2) 



-{AijCiCj + 2BijXiej + DijXiXj) 



where the Einstein convention on repeated indices has been used, and the first 
and third terms are the stretching and bending energy densities while the second 
term characterizes the energy associated with the coupling between bending and 



stretching. Since it is often useful to know the constitutive equations (2.1) in their 
semi- inverted form, we write them as: 



a 


/3 




S 




e 

M 



t We use Reissner's notation for the twisting curvature, hixy 
height of the middle-surface above the xy plane. 



(2.3) 



-2d^w/dxdy, where w is the 
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where: 

a = A-S /3 = -A-^B, ^ = D-BA-^B. (2.4) 



so the energy density (2.2) can be expressed in terms of N and x ^s: 



^ = l{(^^JN^NJ + S^jXzXj) . (2.5) 

In absence of body forces or external loads, the conditions for kinematical compat- 
ibility and mechanical equilibrium translate respectively into the following partial 
differential equations (Calladine 1993): 



dy'^ dx dy dx'^ 



-SK, (2.6a) 



d'^M d'^M d'^M 

-d^ + ^ + ^ = ^^^^ + ''^^^^^ + ""^^y ' (^-^^^^ 

where 5K = K — Kq is the difference between the Gaussian curvature of the 
deformed plate K = z.^ and its spontaneous Gaussian curvature Kq = 

CxCy — Cxy The radius of curvature of the strip is assumed to be much smaller than 
the typical length scale of the deformations (i.e. shallow shell approximation). This 
allows to simplify the differential structure of the theory by transforming covariant 



derivatives into standard partial derivatives. Eqs. |2.6| are then the classic Foppl-von 
Karman plates equations expressed in terms of stress and curvature. Since the strip 
is assumed to be infinitely long in the y direction, all physical quantities should 
be invariant with respect to translations along y. Mechanical equilibrium requires 
V- N = 0, which together with the condition of translational invariance along y 
yields: 

dN, ^ dN,y ^ ^ 

dx dx 

In absence of applied forces along the edges this implies = N^y = 0. In addition, 
the Mainardi-Codazzi compatibility equations Vi^jk = ^k^ij yield: 

dy dx ' dx dy ^ 

that, together with the condition of translational invariance along y, imply that Ky 
and hixy are both constant. Thus, equations (2.6) simplify to: 



^ = -«, (2.7a) 



Using the semi-inverted constitutive equations ( |2.3[ ), one can easily rewrite ( |2.7| ) in 
the form of a single differential equation for the transverse moment M^. To do this 
we start by writing: 

Cy = a22Ny + f32iXx + f^22Xy + f^26Xxy , (2.8a) 



Mx = -p2lNy + hlXx + h2Xy + hQXxy • 



(2.8b) 
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Then, inverting (2.83) and using ( |2.7| 3) we get: 

1321 d^M^ 



— '-a; 



1 

dx'^ 5ii \ dx'^ ' tZy dx^ 



H2Xy - OiQXxy , 



1 fd^M, , p2id^M, 



(2.9a) 
(2.9b) 



Combining these two equations with (2.7) we finahy obtain: 



= {Sll{C^Cy - + Kly) + Ky[Sl2{Ky - Cy) + 5ie{K^y 



hlCx]} ' 



(2.10) 



The boundary conditions associated with equations (2.7) for a strip with no apphed 
forces and torques at the lateral edges have been discussed by Reissner (1993) and 
translate into the requirement: 



dx 







(2.11) 



Together, (2.10) and (2.11 ) complete the formulation of the boundary value problem 



for the behavior of a long anisotropic strip with spontaneous curvature and twist. 

(b) Anisotropy classes 

Before we go further, it is useful to review the different form of elastic anisotropy, 
following the theory of laminates (Jones 1999). Various forms of matrix coupling 
the in-plane strain to the out-of-plane curvature B can be obtained by controlling 
the relative orientation of the laminae forming the layers of a composite plate. 
Thus in a cross-ply fiber reinforced composite, in which the fibers in each layer are 
alternatively oriented at 0° and 90° with respect to the y axis of the strip, one has: 



(2.12) 



with Bii = E±_h?{E\\/Ei_ - 1)/4P, where E\\ and Ei_ are the Young moduli in 
the parallel and transverse directions of the fiber and P is the number of layers of 
the laminate. Antisymmetric angle-ply laminates, on the other hand, have laminae 
oriented as same angle 6 with respect to the laminate coordinate axis on one side of 
the middle-surface and the corresponding equal thickness laminae oriented at — 6> 
on the other side at the same distance form the middle-surface. In this case: 














-^11 





















B16 












B26 






(2.13) 
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Bl2 







B22 











^66 



General asymmetric laminates obtained by the asymmetric stacking of isotropic 
layers with different material properties about the middle-surface, have finally: 



(2.14) 



In the next section we start with a focus on the case of orthotropic strips for which 
Bij = 0, while the case of angle-ply laminates is discussed in Section [4] in the context 
of tristable strips. 

(c) Equilibrium conformations 
In the case of naturally flat strips with no coupling between stretching and 



bending (i.e. Bij = and 5ij = i^ij), equation (2.10) reduces to the equation given 
by Reissner (1993) for the case of a inhomogeneous anisotropic strip: 



OL22 ^ . + = ^ 



Now, letting: 



P2ll^y 

(y-22hi + P21 



D 



11 



^2 



Ot22hl + P2I 



and: 



9^3 



S {hl{CxCy - + H.ly) + t^y[8i2{l^y - Cy) + (^i6(/^a:y " ^^xy) " ^^llC^:]} 

(^22^1 + P2I 



the solution of equation (2.10) which satisfies the homogeneous boundary conditions 



(2.11 ) is given by: 

= — ^^ii(^ — ^0) + ^^iiC'i cosh kix cos k2X + (^iiC2 sinh kix sin k2X , (2.15) 



where: 



{CxCy Cxy) ~^ ^xy 



and: 



with: 



Ci = {hi - hio) 
C2 = {k.- k.o) 

ki-- 



^x ~r ^ yCy Kjy ) \ r- \^xy ^xy ) • 

On On 



ki cosh aki sin ak2 + k2 sinh aki cos ak2 
ki sin ak2 cos ak2 + /c2 sinh aki cosh aki 

k2 cosh aki sin ak2 — ki sinh aki cos a/c2 
/ci sin ak2 cos a/c2 + /c2 sinh aki cosh a/ci 

k2 = 



(2.16a) 
(2.16b) 

(2.17a) 
(2.17b) 
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The longitudinal stress Ny and the transverse curvature hvQi^ can be readily calculated 
by replacing in equations ( |2.7| 3) and ( |2.9^ ). This yields: 

= ^iiC2 cosh kix COS — SiiC[ sinh kix sin , (2.18) 

with: 

^/ _ (^22^11)^^1 P21C2 ^, _ (0^22^11 )^g2 - hiCi 
and finally: 

= ^0 + {C\ + /321C2) cosh /cix cos k2X + (C2 — hiC'i) sinh /cix sin /c2X . (2.19) 



Eq. (2.19) is equivalent to an expression given by Galletly & Guest (2004b) for the 
case of a bistable composite slit tube. We will now examine the effect of spontaneous 
curvature and elastic anisotropy on the energy landscape of elastic strips. 



3. Multistable configurations of orthotropic strips 



As a starting point in our analysis we consider the case of orthotropic strips. For 
this class of materials there is no coupling between stretching and bending, hence 



— and the A and D matrices can be expressed in the form: 

r 1 77 1 

t^E^ 



A = tE^ 



1 V 
p 



D = 



12 



1 V 
p 



where £^0 = E / {1 — j ^) and, following Seffen (2007), we have called: 



E:r=E, 



Ey = pE, G = pE, 



^yx 



^xy 



13' 



where Ex and Ey are the Young's moduli in the x and y direction, G is the shear 
modulus and u^y and Uy^ the Poisson ratios. The conditions for the matrices A and 
D to be positive definite, translate into the requirement P > u'^ and p > (Vidoli & 
Maurini, 2009). The transverse curvature (2.19) and the longitudinal stress ( 2.18| ) 
then simplify to yield: 



^0 + Ci cosh kx cos kx + C2 sinh kx sin kx , 



Ny = {tpED) 2 {C2 cosh kx cos kx — Ci sinh kx sin kx) , 



(3.1a) 
(3.1b) 



where D = t^E/l2(l - v^l^), ki = k2 = k and: 

The stretching energy density Us and the bending energy density u^ are obtained 
by substituting (3.1) into (2.5) and yields: 

1 



1 



D 



2tpE y 



+ PXy + "^^XxXy + P 1 



1 



A-xy 



(3.2a) 
(3.2b) 
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per unit length: 



Integrating (3.2) along the width of the strips we obtain the total stretching energy 

(3.3) 

cosh 2ka — cos 2ka ^ sinh 2ka sin 2ka 



where: 



sinh 2ka + sin 2ka ' 
and the bending energy 



(sinh 2ka + sin 



Ut = aD 



where: 



7/2 



Cy)'^ + 2u{K.O - Cx){l^y - Cy) 



1 



Cxy) + ^(«; - Ko)* + (t - to) * 



(3.4) 



(3.5) 



Before proceeding with the analysis, it is useful to introduce a set of dimensionless 
quantities defined by: 



X 

a 



t 



acj 



{i = X, y, xy). 



so that the dimensionless stresses and energies per unit length are given by: 



TV, 
aE ' 



taE ' 



taE ' 



U = Us^U,. 



U is thus the total elastic energy of a strip that accounts for both stretching and 
bending. In the following sections we will analyze the "energy landscape" embod- 



ied in U and given by (3.3) and (3.4) in a variety of scenarios. We will start by 
considering configurations of constant transverse curvature kx and then move on 
the more general case in which kx is allowed to vary across the strip. In section ([c| 
we will examine the limit of vanishing bending stiffness and show how tristability 
arises in this setting for a broad range of spontaneous curvatures. 



(a) Uniform curvature configurations 
To understand the energy landscape of the strip, we analyze the total elastic 



energy U given by (3.3) and (3.4). We first focus on configurations having constant 



transverse curvature kx'- if ky is nonzero, from equation (2.17) and (3.1 i) it follows 
that Ci = C2 = and k = k^. This solution has zero longitudinal stress and, as a 
consequence, is an isometry of the base configuration: a deformation that preserves 
the local metric of the shell so that e = 0. The latter statement can be also verified 
by calculating the Gaussian curvature of the deformed strip. Since kx = k = ko^ 



equation (2.16i) yields: 



^xy) ~^ ^xi 



(3.6) 
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which imphes that kx^y — k^y = CxCy — c^y as required by Gauss' theorema egregium 
(Kreyszig 1991). The corresponding total elastic energy per unit length U is given 
by: 



(3.7) 



The condition k = ko used to obtain the isometric solution (3.6) describes, in the 



plane (ky^k^y) a continuous set of configurations of zero stretching energy. The 
explicit form of this geometrical locus is given by: 



-2 I -2 



{Cx + ^^y)^y — ^xy 



(3.8) 



that corresponds to an ellipse with focii at : 

1 



2iy 



{Cx ~^ ^^y) •) f^xy, — 5 



and semi- axes a„ 



cly. The 



and axy = l/i with = ^(c^, - i^Cy) 
bending energy, however, suppresses this soft mode by selecting a minimum in 
this class of isometric deformations: the base configuration having ky = Cy and 



'^xy 



-xy 



A special but not necessarily isometric solution with constant curvature can be 

(3.9) 



found by taking the limit ky ^ ^ oi equation (3.1 1) which gives: 



(3.10) 



The corresponding longitudinal stress is given then by: 

^v = \l^ (3*' - Mc^Cy - c%) + kly] 
using which we can calculate the total elastic energy as: 

;^2 X , ^2 12 , ^ 
■^xy) ^ '^xyl ^ -j^2 

This class of deformations belongs to the set of isometrics only when the first term 



} 



(3.11) 



in the left-hand side of equation (3.11) vanishes, namely when: 

,0.2 _ - 



""xy 



^xy 



Co- Co I 



-Ko. 



(3.12) 



Clearly this is possible only when the initial Gaussian curvature Kq is negative or 
zero. In the latter case the strip has only one minimum corresponding to the base 
configuration (see right of Figure [5]). 

However, a secondary equilibrium configuration of constant curvature is ob- 



tained if Cx = Cy = and Cxy ^ 0, when Eq. (3.11) reads: 



u = p 



{45 ^^'y 



(3.13) 



Then, in addition to the base configuration kxy = Cxy^ the energy has a secondary 
minimum corresponding to the twisting curvature: 



^xy 




15 



(3.14) 
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Figure 2. Contour plot of the total elastic energy per unit length U from (3.3) and (3.4) 



in the plane [ky^ kxy)- The labels P and S indicate the primary and secondary minima. 
On the left the case Cx — Cy — and Cxy — 2. The strip has a secondary minimum at 
kxy ~ —Cxy. On the right, Cx — Cxy — and Cy — 2. The strip has a single equilibrium 
configuration. 

The other relevant parameters are /3 = 10, p=l, i^ = l/3, f = 10~^. 

which is approximately kxy ~ —c-xy (see Figure [2j left). The stability of the sec- 
ondary minimum depends on p, /3 as well as the Poisson ratio v. To assess this, we 
calculate the components of the Hessian matrix: 



dki p\i 



6 45 



16 {v"^ - 13) 
945 £2 



c2 
xy 



dkydkxy 



= 0. 



dkly 



4 P 



where kxy is given by Eq. (3.14). In order for the secondary minimum to be stable, 



the eigenvalues of the Hessian matrix must be positive. Upon expanding d'^U /dky 
at the second order in t, this condition reduces to the following inequality: 

/2 ^ 



10 

y 



J 



> 0. 



This implies the following requirment on p: 



P< 



10/3-3t/2 
7 



(3.15) 



In summary, strips having zero spontaneous principal curvatures Cx = Cy = and 
non-zero spontaneous twist Cxy^ admit a secondary equilibrium configuration where 
the principal curvatures vanish, i.e. kx = ky = while the twisting curvature 
is reversed kxy ^ —Cxy The secondary minimum is stable for low enough shear 
modulus, a scenario previously discussed by Seffen (2007) for the case of an elliptical 
plate endowed with pure twist. The most important difference between the two cases 
is that here the strip has zero bending moment along the lateral edges unlike the 
plate investigated by Seffen which can support a finite bending moment at the 
boundary. As a consequence the transverse curvature kx^ which is a free parameter 
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Figure 3. Contour plot of the total elastic energy per unit length U from (3.3) and (3.4) in 
the plane {ky, kxy). The labels P and S indicate the primary and secondary minima. On 
the left the case = Cy = 2 and Cxy = 0. In addition to the base configuration, the strip 



has a second equilibrium at ky 



xl/2 



. On the right, Cx 



-2, Cy — 2 and Cxy — 0. 



In this case the strip has only one equilibrium configuration. Note that the region ky <0 
that in the previous example was separated form the region ky > ^ hy an energy barrier 



of height Uo - is now connected by two "passes" located at kxy = ±(c 
other relevant parameters are /3=10, p=l, zy— 1/3, t = 10~^. 



xl/2 



The 



in uniform curvature models, is here adjusted in order to maintain the boundary 
moment-free and this severely restricts the space of possible uniform curvature 
configurations. 

The previous situation is an example of bistability in a strip with a negative 
spontaneous Gaussian curvature Kq = — c^^, where the principal curvatures Cx and 
Cy vanish identically while the surface is purely twisted. However when c^Cy ^ 
and Kq < 0, the strip also exhibits a secondary equilibrium configuration, albeit 
with a variable transverse curvature kx discussed in the next section. 



{b) Non-uniform curvature configurations 

When the constraint of constant curvature is lifted, new minima appear in the 
energy landscape. Figure [S] shows a contour plot of the total elastic energy in the 
plane {kx^ kxy) for the two cases CxCy > (left), CxCy < (right) and Cxy = 
in both case. For positive spontaneous Gaussian curvatures the energy minima are 
those described in the previous section when kxy = 0. The transverse curvature kx 
associated with the secondary minimum dit ky < is mostly constant across the 
width of the strip and equal to 

{CxCy — C^y) + k'^y 

f^o = ^ , 

hZy 

with exception for a boundary layer (Figure [4|. This phenomenon, first noted by 
Lamb (1891), is due to the rapid build-up of the bending moment Mx from zero at 
the edges to a non-zero value that inevitably develops in any configuration other 
than the base state (see Mansfield 1989 for a detailed explanation). An estimate 
of the size of the boundary layer is given by 1/k ^ {t/\i^y\)^] however it vanishes 
in a strip of lenticular cross section whose thickness (hence the bending moment) 
smoothly tapers and vanishes at the edges. 
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Figure 4. (Left) The dimensionless transverse curvature kx associated with the secondary 
equiUbrium configuration of a strip with spontaneous positive Gaussian curvature as a 
function of the dimensionless distance x from the centerUne. The relevant parameters are 
Cx = 2, = 1, Cxy = 0, /3 = 10, p=l and i — 10~^. The transverse curvature is constant 
and equal to hiq with exception for a boundary layer. In the inset a schematic representation 
of the primary (P) and secondary [S) equilibria. (Right) Monostability/bistability (M/B) 
phase diagram in the plane (/3, Cxy) for a strip of negative Gaussian curvature and different 
values of the Poisson ratio The spontaneous principal curvatures are = 2 and Cy — 1. 
The other relevant parameters are p — 1 and i — 10~^. 




Figure 5. The elastic energy given by (3.3) and (3.4) as function of ky and 



for 



strip of spontaneous positive Gaussian curvature (on the left, Cx — Cy — 1) and a strip 
of spontaneous negative Gaussian curvature (on right, Cx — —Cy — 0). In the former, 
the ky < ^ portion of the energy landscape is separated from the ky > ^ region by a 
large energy barrier of height (j{ky — 0) ^ Kq . This implies the existence of a secondary 
minimum (labelled as S). In the case of strips of negative spontaneous Gaussian curvature, 
on the other hand, the energy barrier has two passes at /^^ = and kxy — ±(— ^o)"^^^ 
and the secondary minimum is replaced by a saddle point. The flat plateau in the energy 
barriers is a graphical artifact due to the limited range of U shown. 



On on the other hand, for negative Gaussian curvature CxCy < and Cxy = 0, 
the secondary minimum becomes a saddle point and the elastic energy has only 
one minimum corresponding to the base configuration. It is interesting to notice, in 
this latter case, that the regions of positive and negative ky of the energy landscape 
are separated by a barrier of height U{i^y = 0) ~ Kq that is large everywhere with 
exception for two "passes" at kxy = ±(— c^jC^y)^/^ (see Figure p|. The existence of 
these passes along the energy barrier at = is the reason for the non-existence of 
a secondary equilibrium configuration in strips with spontaneous negative Gaussian 
curvature. As explained in the previous section and summarized in equation (3.12), 
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the stretching energy associated with the barrier at = scales hke Us ~ {Kq + 
kly)'^ and, in presence of spontaneous negative Gaussian curvature, it can be reheved 

through twist. This leads to the formation of passes at kxy = ±(— i^o)^^^ shown 
in Figure [3] and |5] which allow the transverse curvature ky to be switched from 
positive to negative (and vice versa) isometrically. Positive spontaneous Gaussian 
curvatures, on the other hand, cannot be accommodated through twist and thus 
there are no passages across the energy barrier at ky = 0; consequently there is a 
secondary minimum. 



Figure 6. Sequence of deformations of a strip of spontaneous positive (top) and nega- 
tive (bottom) Gaussian curvature associated with curvature reversal. While the former 
sequence requires the strip to become temporarily flat and thus produces large strains, 
the latter can be performed isometrically at zero stretching energy cost. The existence of 
this low energy deformation in strips of spontaneous negative curvature is responsible for 
the formation of the passes shown in Figure [5] hence the lack of a secondary minimum. 



To get an intuitive sense of these passes in the energy landscape, imagine a strip 
of positive spontaneous Gaussian curvature (upper left frame of Figure |6|. To evert 
the strip, one choice is to flatten the strip and then reverse its curvature. Since the 
strip is naturally curved, the intermediate flat conformation will suffer a large strain 
and this corresponds to the large energy barrier shown in Figure |5j On the other 
hand, if the strip is naturally shaped like a saddle (bottom left frame of Figure [6|, 
turning the strip inside-out can be achieved by pulling apart two opposite corners 
until the strip is fully extended and maximally twisted, and then bend it into the 
new configuration. The latter deformation can be performed without stretching the 
surface and thus with no energy cost other than bending, in contrast with the case 
of a strip with positive Kq. Similar deformation pathways have been analyzed by 
Fernandes et al. (2010) as possible strategies to achieve shape control of bistable 
composite plates using embedded actuators. 

However, when strips with negative spontaneous Gaussian curvature < 
have a non-zero spontaneous twist Cxy^ they admit bistable configurations. In addi- 
tion to the configuration with constant curvature described in the previous section 
in the special case Cx = Cy = {) and Cxy 7^ 0, another configuration can be obtained 
by setting CxCy 7^ with i^o < 0. Figure [t] shows a contour plot of the total elastic 
energy (3.3) and ( |3.4[ ) of a strip with c^; = 2, = 1 and Cxy = 2, thus Kq = —2. 
The plot on the left (in which (3 = 10) shows a secondary minimum in the proxim- 
ity of the point {ky^ kxy) = {cy, —Cxy)- When anisotropy parameter /3 is decreased. 
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Figure 7. Contour plot of the total elastic energy per unit length U from (3.3) and (3.4) 
in the plane {k,y, kxy). The labels P and S indicate the primary and secondary minima. 
In both plots the spontaneous curvatures are set Cx — 2^ Cy — 1 and Cxy — 2. On the left, 
the case /3 = 10, in addition to the base configuration, the strip as a secondary minimum 



in proximity of the point {ky 



Such a minimum disappear when ^ — 1 



(right). The other relevant parameters are />=1, z/=l/3, t=10 



however, the secondary minimum disappears and the energy landscape is again 
characterized by a unique minimum corresponding to the base configuration. 

Such a configuration has energy tj ^ {p/j3)c^y^ and has a minimum for large 
values of P/p. One might expect that the value of P/p necessary to have a secondary 
minimum of this kind would increase monotonically with Cxy] however this is in- 
correct. In Figure [4] we show the "phase diagram" of the monostability/bistability 
region for a strip with spontaneous curvatures Cx = 2 and c^y = 1 in the plane 
Cxy). Upon fixing the spontaneous twist Cxy, the bistability region can be ac- 
cessed by increasing as expected. However, for fixed the bistable region is 
accessed by increasing the spontaneous twist. The origin of this behavior is related 
with the height of the energy barrier ai ky = separating the secondary minimum 
from the basin of attraction of the base configuration at /^^ < 0. From equation 
(3.11) we see that the height of such a barrier U ~ PKq. For fixed decreasing 
Cxy (or equivalently \Ko\) will lower the height of the barrier until the basin of 
attraction associated with the portion where ky < of the base configuration will 
merge with that of the secondary minimum. For values of Cxy smaller than this crit- 
ical value, the energy landscape is then characterized by a unique minimum (i.e. 
the base configuration) surrounded by a large C-shaped basin of attraction that 
starts from the base configuration and extends in the region ky < penetrating 
the barrier at = through the passes located at kxy — ±(—^0)^^^ (see right 
of Figure [t]). Once again, the presence of two regions across the barrier at /^^y = 
where the stretching energy drops, plays a crucial role in assuring the existence of 
bistable configurations. 

In summary, orthotropic strips exhibit various form of bistability depending to 
their spontaneous curvature and elastic moduli. The simplest bistable configura- 
tions have constant curvature throughout the width of the strip and correspond to 
isometrics of the base configuration. Non-isometric secondary equilibria have, on 
the other hand, variable transverse curvature kx\ this is approximately constant 
and equal to /^o = (^0 + ^xy)l^y bulk of the strip, but abruptly jumps 

to k = Cx -\- i^{cy — kx) near the edges to relieve the bending moment that forms 
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Figure 8. Phase diagram of the monostabiUty (M, in white) and bistabihty region (B, 
shaded) in the plane (cx, Cy) for strips of various spontaneous twist Cxy. Different shades of 
green correspond to positve (hght) and negative (dark) spontaneous Gaussian curvature. 
For zero spontaneous twist, bistable configuration exist only in the region of positive 
Gaussian curvature. For Cxy / 0, on the other hand, bist ability also occurs in a range of 
negative Gaussian curvature. The relevant parameters are ^5 = 10, y9=l, i/=l/3 and 
t = 10"^ 



in the bulk when the strip is in a configuration different from the base one. For 
strips with positive spontaneous Gaussian curvature, in particular, the secondary 
equilibrium configuration has transverse curvature of opposite sign relative to the 
base configuration. Strips with negative spontaneous Gaussian curvature, on the 
other hand, do not possess a secondary equilibrium configuration unless they have 
sufficient spontaneous twist Cxy The existence of bistability in this case also relies 
on the value of the ratio j3 / p between the Young's modulus along the longitudinal 
direction and the shear modulus. Figure [8] shows a phase-diagram of the monosta- 
bility /bistability regions in in the plane (c^^, Cy) for strips of various spontaneous 

twist Cxy. 



(c) Tristahle configurations 

As we mentioned in the introduction of Section [3j the requirement for the ma- 
trices A and D to be positive definite translates into the conditions P > u'^ and 
p > 0. In this section we show that when the limit P ^ u'^ is approached, a new 
stable minimum appears in the energy landscape of the configurations allowing 
for tristability. The presence of three minima in the reduced quasi-one-dimensional 
geometry of elastic strips complements the appearance of such a landscape for cor- 
rugated strips, where tristability arises as the consequence of the coupling between 
internal prestresses created while imprinting the corrugation and non-linear geo- 
metrical changes during deformation (Norman et al. 2008). 



Letting /3 = z/^ -h e and expanding equations (|3.3|) and (3.4) in powers of e gives: 



12 



{ky - Cyf ^ {kly - cly) [>+o(e). (3.16) 



Then, with kxy = Cxy^ it is easy to prove that dtj /dky = when ky = Cy and also 
when: 

^^ = \[^y^l 3c, ± ^ {cx -lycy)^- 30 P | . (3.17) 
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Figure 9. Contour plot of the total elastic energy per unit length U from (3.3) and (3.4) 
in the plane {ky^ kxy) of a strip with /3 = + 10~^ and spontaneous curvatures = 1, 
Cy — 2 and Cxy — 0. The primary, secondary and tertiary minima are indicated with the 
letters P, S and T respectively. The plot on the right shows a magnification of the region 
ky G [1.8, 3.2] where the primary and tertiary minimum are located. The other relevant 
parameters are p=l, z^=l/3, t = 5x 10~^. 



Thus, in addition to the base configuration and the secondary minimum at ky < 
discussed in section (|b|, the energy has a third minimum at while 
ky is the location of the maximum separating the third minimum from the base 
configuration. For relatively thin shells, £ <C 1 so that one finds simply k'^ ^ Cx/i' 
and ky {cy ^ Cx/h') /2. The energy of the third minimum is given approximately 
by: 

U{k^y)^Y2^^'{cy-"f)' (3.18) 

Figure [o] shows a contour plot of the total elastic energy per unit length U from 
^% and for = 1, = 2, c^y = and /3 = + 10"^ with v = 1/3 and 
p = 1. The primary, secondary and tertiary minima are indicated with the letters 
P, 5* and T respectively. 

In the limit (3 ^ u'^, tristability has been noted by Vidoli & Maurini (2009) in the 
case of plates with free boundaries and uniform curvature, although the physical 
origins of this phenomenon were not discussed. In our simple setting, consistent 
with the arguments in section ([a|), there is a continuous set of configurations of 
zero stretching energy corresponding to the isometrics of the base state. In the 
plane {ky^ kxy) this set describes an ellipse whose intersections with the Hvy axis 
(i.e. the untwisted configurations) are given by ky = Cy and ky = Cx/i^. The small 
but finite bending stiffness along the principal directions prevents this from being 
a true zero energy mode and instead leads to the existence of a soft mode by 
raising the energy of the minimum at ky = Cx/i^; for {3 ^ v'^ (i.e. D22 ^ ^11) 
the minimum is completely suppressed. Taking the limit {3 ^ v'^ corresponds to 
lowering the bending stiffness so that the energies of the two isometric configurations 
ky = Cy and ky = Cx/i' are comparable again. This mechanism is completely general 
and independent of the specific geometry of the problem; any plate with shallow 
spontaneous curvature is amenable to a set of isometrics and multistability can be 
in principle obtained by setting det(D) ~ 0, which in the case of orthotropic plates 
corresponds in fact to {3 v'^ . 
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Figure 10. (Left) Total elastic energy as a function of ky for a strip with — — 10~^ and 
spontaneous curvatures = 1, = 2 and Cxy = 0. The inset shows a magnification the 
boxed region ky G [1.8, 3.3]. The height of the barrier separating the base configuration 
from the tertiary minimum located at ~ Cx/z^ = 3 is four order of magnitude smaller 
than the barrier at /^^ = 0. (Right) phase-diagram in the plane (c^, e = /3 — z^^) of the 
bistable (M, in white) and tristable (T, shaded) regions of a strip with zero spontaneous 
twist and various Cy. The other relevant parameter are p — 1 and i — 10~^. 

An interesting feature of this tristable energy landscape relies on the fact that 



the energy barrier separating the minima at ky 



and ky = c^jv is generally 



much smaller than that located at /^^ = separating the base configuration with 
the minimum at ky < 0. More precisely: 

1 



720 



while U(fS) 



(1/45) {vKqY . Assuming kxy 
of these two energy barriers can be approximated as: 

2 

2 



(c, - ucyf (3.19) 
xy = the ratio between the height 



[/(O) 



1 

16 



1 



sinh - loff 
2 ^ 



(3.20) 



Assuming Cx < Cy^ this ratio is O(10~^). For Cx > Cy the previous relation is no 
longer valid as the bending energy, that was neglected to derive ( |3.20 ) becomes 
relevant; however, the order of magnitude of the ratio U{ky)/U{0) remains always 
of order 0(10-^). 

The stability of the tertiary minimum at ^ Cxjv depends, in general, on 
the value of the Poisson ratio and the spontaneous curvatures Cx and Cy. Figure 
10 shows the bistability/tristability regions as a function of e = /3 — and Cx for 
various values of Cy. Increasing the transverse spontaneous Cx allows tristability to 
occur at larger values of e, thus further away from the limiting condition j3 = v'^ . 



4. Strips with coupling between stretching and bending 

We now consider generally anisotropic strips, in which stretching and bending defor- 
mations are coupled through the B matrix, so that the total elastic energy density 
reads: 



+ ^ll{l^x - CxY + '^^12{t^x - Cx){l^y - Cy) 

+ h2{l^y - CyY + ^m{l^xy 



^xy 



f] (4.1) 
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where the elastic moduh a22 and Sij are given by (2.4) and Ny and have the 



general form given in equations (2.18) and (2.19). Because of the reduced dimen- 



sionality of this class of elastic objects, the coupling between stretching and bending 
is not expected to produce dramatic differences when compared to the case of or- 
thotropic strips discussed in the previous section and indeed, for the majority of 
the choice of the elastic moduli a22 and Sij^ the behavior is qualitatively identical 
to that already described. 

Under some circumstance, however, the finer control of the material proper- 
ties characteristic of anisotropic laminates might be taken advantage of to enhance 
the region of tristability. To understand this possibility, we first note that the re- 
quirement for the matrices a and 6 to be positive-definite, translates here into the 
condition 0^22 > 0, ^66 > and S11S22 > ^12- If the latter condition, in particular, is 
only weakly satisfied, the strip will be tristable for some value of the spontaneous 
curvatures Cx and Cy. To illustrate this point let us consider the case of an antisym- 
metric angle-ply laminates. The condition ^11^22 > ^12 translate into the following 
inequality: 



(^16 ~ ^66^1l)(526 ~ ^66^22) > {B1QB26 — ^66^12)^ . 



Regardless the specific values of Bij and Dij^ the left and right hand sides of this 
inequality will be comparable in any material where ^ AqqD. In other words the 
limit det(^) 0, that for orthotropic strips can be approached by taking j3 ^ v'^ ^ 
can here be obtained by increasing the strength of the coupling between stretching 
and bending. Thus we might expect angle-ply laminated strips with strong coupling 
between bending and stretching to be tristable for a broad range of spontaneous 
curvatures. 



5. Conclusions 

Despite its relatively long history, the theory of elastic mulitstability is still a field 
rich in challenging problems and open questions. Most of the theoretical work has 
been limited to the case of uniform Gaussian curvature deformations. Here, we 
lift the assumption of uniform curvature variations and inextensibility, and find a 
variety of morphing and multistable scenarios in shells that have both anisotropy 
and spontaneous curvature. Our results complement and extend those previously 
reported by Galletly & Guest (2004a, 2004b), Guest & Pellegrino (2006), Seffen 
(2007) and Vidoli & Maurini (2009). In particular we emphasize the fundamental 
role of the intrinsic geometry of the base configuration embodied in the spontaneous 
Gaussian curvature of the strip. Both strips of positive and negative spontaneous 
Gaussian curvature admits bistable configurations, but while strips of positive spon- 
taneous curvature are always bistable, regardless the values of the elastic moduli, 
the occurrence of bistability in strips of spontaneous negative curvature depends 
crucially on the presence of spontaneous twist as well as the relative stiffness of the 
strip under tensile and shear deformations. This fundamental difference stems from 
the fact that the principal curvatures of a saddle-like strips can be isometrically 
switched from positive to negative and vice-versa. This leads to the formation of 
low energy pathways across the energy barrier corresponding to the fiat configura- 
tion. 



Multistability of free spontaneously- curved anisotropic strips 



19 



Furthermore, we find that tristabihty occurs in strip-like plates when det(J) ^ 0, 
where S = Y> — BA~^B is the effective bending stiffness matrix. In orthotropic ma- 
terials, this is corresponds to the limit /3 ^ u'^ where P = Ey/E^ discussed by 
Vidoli & Maurini (2009) for the case of plates with uniform curvature. A special 
feature of tristable strips is associated with the height of the barrier separating 
the two equilibrium configurations of like-sign curvature. The height of this bar- 
rier is, in general, much smaller than that associated with the fiat configuration. 
This might have interesting implications for nanoscale materials. Typical strip-like 
biopolymers, such as proteins or cytoskeletal filaments, have widths of the order of 
few nanometers, lengths that ranges from 1 to 10^ nm and elastic moduli of the or- 
der of mega-Pascal. Assuming c^^, Cy and Cxy of order one and at ^ 1 nm^ one finds 
that U{Ky) ~ 3 • 10~^ /c^T/nm. Thus even a micron-long strip-like polymer with 
non-zero spontaneous curvature would be able to fluctuate between minima under 
the sole effects of thermal excitations. Extending our analysis of multistability to 
include the effects of thermal fluctuations is likely to yield insights on the dynamical 
behavior of biopolymer assemblies and is but one extension of our present analysis. 
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